Quasi-normal modes of charged, dilaton black holes 



Valeria Ferrari^ 1 )*, Massimo Pauri^ 2 ^ and Federico Piazza^ 3 ^ 
^Dipartimento di Fisica, Universitd di Roma, "La Sapienza" and 
I. N.F.N. Sezione di Romal, Piazzale Aldo Moro 5, Roma, Italy 
Dipartimento di Fisica, Universitd di Parma and 
I. N.F.N. Gruppo Collegato di Parma, Viale delle Scienze, 43100, Parma, Italy 
Dipartimento di Fisica, Universitd di Milano and I. N.F.N. Sezione di Milano, Via Celoria 16, 20133 Milano, Italy 

(Received 30 May 2000) 



In this paper we study the perturbations of the charged, dilaton black hole, described by the 
solution of the low energy limit of the superstring action found by Garfinkle, Horowitz and Stro- 
minger. We compute the complex frequencies of the quasi-normal modes of this black hole, and 
compare the results with those obtained for a Reissner-Nordstrom and a Schwarzschild black hole. 
The most remarkable feature which emerges from this study is that the presence of the dilaton 
breaks the isospectrality of axial and polar perturbations, which characterizes both Schwarzschild 
and Reissner-Nordstrom black holes. 

PACS numbers: 04.30.Db, 04.70.Bw, 04.40. Nr 



I. INTRODUCTION 

The study of the effects produced by the coupling between gravity and scalar fields is a rather hot issue, since a 
definite prediction of superstring theory is just the existence of a scalar field, namely the dilaton. For instance, this 
coupling could have played an important role during the early phases of the life of the Universe. Indeed, it has been 
shown that a possible cosmological implication of the dilaton is the so-called "pre-big-bang" scenario Q , characterized 
by a dilaton-driven inflationary period, that should have left its peculiar fingerprints on the cosmic gravitational wave 
background ||. On the experimental side, recent investigations on the interaction of scalar waves with gravitational 
detectors have shown that a scalar component of the gravitational radiation should excite the monopole mode of a 
resonant spherical antenna Q and also give rise to specific correlations between the signals revealed by a resonant 
sphere and an interferometer ||. 

Among all possible astrophysical sources of gravitational waves, black holes should have the most typical and in 
principle recognizable frequency spectrum. According to the Einstein theory of gravity, a black hole which dynamically 
interacts with its surroundings emits wave bursts having shape and intensity which initially depend on the features of 
the external perturbations. As a late-time effect, however, it generates gravitational wave trains having characteristic 
frequencies which are independent of the initial perturbations, the so-called quasi normal modes ||. And we know 
from the theory that these frequencies can only depend on the three parameters characterizing a black hole, namely 
the mass, the charge and the angular momentum. 

The question whether a black hole can support an additional scalar degree of freedom and still generate a regular, 
asymptotically flat, spacetime free of naked singularities or any other kind of pathological behaviour has been widely 
investigated in the literature. Having in mind the "no-hair theorem" program, Bekenstein |7j in 1972 showed that 
black holes cannot support a minimally coupled scalar field, if the solution is required to be asymptotically flat. 
Subsequently, he developed a procedure to generate black hole solutions surrounded by a conformal scalar field 
||, but they proved to be unstable against radial perturbations ||. In 1990 Ferrari and Xanthopoulos flOfl , using a 
Kaluza-Klein approach, derived equations describing a gravitational field coupled with a massive scalar field, obtained 
the conformal structure of the metric, and discussed several kind of possible couplings. Yet, no definite answer was 
given to the question whether a satisfactory black hole can exist surrounded by a massive scalar field. 

The search for regular exact solutions describing black holes endowed with a scalar field got a formidable boost when 
it became clear that, as long as the curvature is small, all vacuum solutions of general relativity are also approximate 



*Email address: Valeria.Ferrari@romal.infn.it 
t Email address: pauri@pr.infn.it 
"'"Email address: piazza@mi.infn.it 



1 



solutions of string theory. As a matter of fact, charged black hole solutions coupled to a scalar field, viz a dilaton in 
this context, can be obtained by the action (we use geometric units throughout) 

S = J d 4 x^g~ [R - 2(V0) 2 + e- 2a *F\ (1) 

which, by variations, gives the equations: 

V^e-^F^) = 0, (2) 

V 2 </> - i e - 2a ^ 2 = 0, (3) 

R»u = 2V^V„0 - 2e- 2 ^F^Ff + ig^e'^F 2 , (4) 

where F 2 = F fl ^F^ lJ is the first Maxwell invariant and a is a non- negative real constant regulating the strength of 
the coupling between the dilaton and the Maxwell field. This theory has a number of interesting limiting cases. The 
limit a = corresponds to the ordinary Einstein-Maxwell theory plus a Klein-Gordon scalar field with zero mass. The 
case a = \/3 is the four-dimensional reduction of Kaluza-Klein theory. Finally, for a = 1 the action (|l|) describes the 
tree-level low energy limit of superstring theory in the so-called Einstein frame. A large class of black hole solutions of 
this theory in an arbitrary number of dimensions has been found by Gibbons and Maeda [ fi"i"l . Their results were next 
specialized to four dimensions by Garfinkle, Horowitz and Strominger [Q (GHS) who studied a solution describing a 
spherically symmetric, charged, dilaton black hole. In the GHS solution the electric charge and the dilaton are not 
independent parameters: when the charge is set equal to zero the dilaton also disappears, and the solution reduces 
to that of a Schwarzschild black hole. In a sense, this is a consequence of the "no-hair theorems" which limits the 
number of free parameters of a black hole to three. 

The equations satisfied by small perturbations of this solution have been derived and studied by Holzhey and 
Wilczck jl3| in the line of the general treatement given by Chandrasekar They first showed that it is possible 
to reduce the perturbed equations to five decoupled wave equations with potential barriers, and, as a by-product of 
their analysis, they argued that the GHS solution is stable under external small perturbations. 

In this paper we take a further step in the study the properties of the coupled emission of electromagnetic, scalar 
and gravitational radiation by such black holes, computing the quasi-normal mode frequencies of the GHS solution. 
Comparing the spectrum of the latter with those of Schwarzschild and Reissner-Nordstrom, we show that the presence 
of the scalar field breaks the relevant feature of the isospectrality of the axial and polar perturbations. 

The paper is organized as follows. In Section 2 we summarize the features of the GHS solution. In Sections 3 and 
4 the equations governing the axial and polar perturbations are discussed separately. We shall not derive the axial 
equations, since there is nothing to add to the derivation made in Ref. J13J . On the other hand, the polar equations 
will be discussed in greater detail. Some misprints appearing in Ref. |13[| are corrected and the explicit expression of 
the matrix whose eigenvalues are the potentials governing polar perturbations is given in the Appendix. In Section 5 
the quasi-normal frequencies of the GHS black hole obtained with an extended WKB approch |L6| are calculated for 
different values of the parameters. Finally, in Section 6 the results obtained are discussed and the main qualitative 
differences between the spectra of the dilaton and the Reissner Nordstrom black holes are pointed out. 

II. THE EXACT SOLUTION 

The exact solution of equations (||)- ([!]), with a = 1, describing the charged dilaton black hole we study in this paper 
is § 



ds' 



1 _ d e -(\- dr 2 -r(r- Q±j [d9 2 + sin 2 6dv 2 ] (5) 



where M and Q e are, respectively, the black hole mass and electric charge. F tr = Q e /r 2 is the only non-vanishing 
component of the electromagnetic tensor, and the scalar field is related to the electric charge by the following equation 

e " - f 1 " I) ■ < 6 » 
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which shows that vanishes at radial infinity. Since the dilaton and the electric charge are coupled through Eq. 
@ , the Reissner Nordstrom solution cannot be obtained as a limiting case of the metric @ . On the other hand, the 
Schwarzschild solution is recovered from Eq. (||) by setting Q e = 0. 

It should be noted that the usual relation between radius and area of the spheres t = const, r = const is obtained 
in terms of the modified radial variable 




and not of r, by the usual relation A(r) — Airr 2 . 

The metric ([3j) appears to be singular in r = 0, r = Q^/M and r = 2M. The surface r = 2M is an event horizon, 
whereas on r — Q 2 /M the curvature scalar 



2 V r ) r(Mr - Qf) ' 

diverges, showing that there is a curvature singularity. In r = Q^/M, the radial coordinate f vanishes and loses 
its meaning for r < Q 2 e /M. Thus, a physical observer who crosses the event horizon terminates its jurney on the 
curvature singularity. 

In order to avoid naked singularities, in what follows we shall assume that 

2M 2 > Q 2 e . 

As shown in Ref. |T^] (to be referred to hereafter as MT), the study of the perturbations of the metric (||), can be 
restricted, without loss of generality, to axisymmetric perturbations only. The appropriate metric in this case is 

ds 2 = e 2u {dt) 2 - e 2 ^{d(t> - q 2 dx 2 - q 3 dx 3 - udt) 2 - e 2 ^{dx 2 ) 2 - e 2 ^ (dx 3 ) 2 , (7) 

where, in the unperturbed state, 

e 2u = e -2 M2 = ^ _ t ^ = ~2 sin 2 Q 

e 2M3 = r 2 , lu = q 2 = q 3 = , (8) 

and (0,1,2,3) stand for (t, ip, r, "!?). We shall assume that, as a consequence of a generic perturbation, the metric 
functions, the electromagnetic quantities and the scalar field will experience small changes with respect to their 
unperturbed values 

V > V + 5v, [12 >U2+S^2: ^ — >ip + Sif), 

(l 3 > [i 3 + Su 3 , u > Suj, q 2 > Sq 2 , 

q 3 — > Sq 3 , — > F^ u + 5F^ Vl (f> ► (f> + 8cj), 

Since the perturbation is assumed to be axisymmetric, all perturbed quantities depend on t, r and 9 only. As in MT, 
it is convenient to project Einstein's and Maxwell's equations onto an orthonormal tetrad frame, and assume that 
all perturbed functions have the time dependence e lCT *. The separation of variables is accomplished by expanding all 
perturbed tensors in tensorial spherical harmonics. These harmonics belong to two different classes, depending on 
the way they behave under the angular transformation 6 — > it — 9 and <p — * n + </>. Those that transform like 
(— are termed axial, those that transform like (— 1)W are termed polar. The perturbed equations split into 
two decoupled sets corresponding to a different parity. 



III. THE AXIAL EQUATIONS 

The axial equations for £ > 2 are obtained by perturbing the {12} and {13}-components of the Einstein equations. 
The separation of variables, can be accomplished by putting 

fV sin 3 % 2)3 (t,r,tf) - <fe, a (t,r,0)] = Qt(r, a)C^' 2 (tf)e^, 
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and 



F ol (t,r,0) sini? = 3B,(r,a)C7 + Y 2 (i9)e lCT *, 

where C;, 1 / 2 ^) are the Gegenbauer polynomials. As shown in Ref. the axial equations can be cast in the 

following form 



(ir 2 



Hie 



B 



Hit 



where 



a 2v 



B 



2 + 2 + ^ 4- 3< ^ e r 2l/ 
r 2 4M 2 f2 









) i 







2/LtQ e -6M 



r* is the tortoise coordinate defined by the equation 



~ 2v dr. 



(9) 



(10) 



(11) 



and [j, 2 = (£ — 1)(£ + 2). The radial functions H\ and i?2 (from now on we omit the index t) are related to the 
perturbed metric and electromagnetic functions by the following relations 



Q(r, a) = rH 2 (r, a) 



re" B(r, a) 



ffi(r, a) 
2[i 



(12) 



The right hand side of Eq. (^) can be obviously diagonalized by a linear r-independent transformation, in the form 

(13) 



where 



4M 2 

Then the system (^) decouples in two wave equations 

d 2 



Z" — C\Hi + C2H2, 
Z$ = C2H1 — C1H2, 



9M 2 + Ql (3 + 4^i 2 ) and C 2 = 2nQ e 



drj 2 



a 1 \ Z? = VfZ? 



(* = 1, 2), 



(14) 



(15) 



where the explicit form of the effective potentials is 



„2// 



v a 

v l,2 



(H 2 + 2)r + ^ + 



'ZQjr | 
4M 2 f 2 ' 



2A/ 



3M± 



9M 2 + Q2( 3 + 4ai 2) 



4M 2 



(16) 



Equation ( p"4| ) shows that when Q e vanishes, £2 = 0; in this case Z% reduces to the gravitational perturbation and V 2 a 
to the Regge- Wheeler potential, whereas Zf reduces to the pure electromagnetic perturbation of a Schwarzschild 
black hole, which is known to be independent of gravitational contributions. If the electromagnetic charge does 
not vanish, the gravitational and electromagnetic perturbations are coupled; a gravitational wave incident on the 
potential barriers induces the emission of electromagnetic radiation and viceversa. On the other hand, the dilaton is 
not dynamically coupled with the axial perturbations. Its effect is that of shaping the effective potentials together 
with the electric field. In a similar manner, the energy density and the pressure of the matter composing a perturbed 
star determine the potential barrier of the axial perturbations without being dynamically coupled to the perturbed 
gravitational field. 

The I = 2 potentials V± and V 2 a are plotted in Fig. ?? versus the rescaled tortoise coordinate r*/M for different 
values of the electric charge. We see that they always tend to zero at radial infinity and at the black hole horizon, 
except when the charge assumes its extremal value Q 2 = 2M 2 . In this case they take the form of a step, reflecting 
all waves whose square frequency is lower than their limiting value on the horizon. 
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We also see that for non extremal black holes the maximum of the barrier moves towards the horizon as the electric 
charge increases. 



n — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — q 




" u I i i i i I i i i i I i i i i u " u I i i i i I i i i i I i i i i u 

-10 10 20 -10 10 20 

r,/M r/M 

FIG. 1. The I = 2 axial potentials V\ (left) and V 2 a (right) are plotted versus the rescaled tortoise coordinate for different 
values of Q = Q e /M. When the charge is set equal to zero, V" reduces to that of the wave equation governing the pure 
electromagnetic perturbations of a Schwarzshild black hole, whereas V 2 a reduces to the Regge- Wheeler potential. 



IV. THE POLAR EQUATIONS 



The separation of variables for polar equations is accomplished by requiring that the perturbed functions have the 
angular dependence deriving from the expansion in tensor spherical harmonics: 



N e (r)P e (9) 



Aat 



5v 

Sfi 2 = L e (r)P e (6)e 
5*P = [T e (r)P e (6) + 2fi- 2 X e (r)P e (6), e cab 6] e 
6[i 3 = [T t (r)P e (6) + 2fi- 2 X e (r)P e (8), e ,e] e l °\ 

-B mt {r)P t {9)e M 



icrt 



6F 02 
SF23 



2Qe 

-L-B me (r)P e (9),ee M , 



te ar 

2Q e 

<f>P e (8)e lrTt 



B 2 3e(r)Pe(6),ee 



iat 



where Pg(Q) are the Legendre polinomials. In the following, we shall omit the index I in all radial functions. Among 
the eight radial functions, N, L, T, X, -B02, ifoj, B23 and the function T can be eliminated from the perturbed 
equations by making use of the relation that follows from the {03}-component of the perturbed Einstein equations 



B 23 ~T-L + 2^ 2 X = 0. 



(17) 



From the polar components of the Maxwell's equations (see MT, Chap. 5 Eqs. (165)-(167) for the Reissner Nordstrom 
case), it is easy to derive a second order equation for B23: 



re 2v (r 2 B 23 ) !r 



+ [a 2 f 2 e- 2v 



- ( M 2 + 2)] B 23 



•5 



2Qt 



(2B 23 - 3L - 2X - N - 2$) = 0. 



The Einstein equations for <5i?o2, SR23, SG22 and SRn give: 



_ + (logfe- y ) 



{B 23 -L-X)-(logr), r L 
2 



(TV - i), r - (log re") r L ~ (log f e" 1 ') r N ~-B 23 + 20 r $ 

r 



X, 



2(logre") tr X tT 



2(logf), r iV iI . + 2(logfe 1/ ), I .(B 2 3 
(^ 2 + 2)e- 2 " 



2f 2 
L-X) 



(N- 



L) + a 2 e~ Av X 



(18) 

(19) 
(20) 
(21) 



-2a 2 e- 4 ^(B 2 3 



N - 2(logf), r (logfe 2l/ ) :rJ L + 2(^ r ) 2 L 
2g 2 e- 2 - 



L-X) 



(22) 



Finally, the perturbed equation for the dilaton is obtained by perturbing eq. (||): 



l(f 2 e 2 "$ ,.),,. + (aV 



(p 2 + 2) , 2Q 2 



e IV {N -ZL-2X + 2B, 



23),r(P,r 



, (r 2 e 2 ^ r ), 



^»2 ^ 2 
2;/ j 



02 



0. 



(23) 



The seventh order linear system composed by equations (^8|)-(p3|) has been shown to be reducible to three Schrodinger- 
like equations by using the following procedure. It is well-known that the order of a system of linear differential 
equations can be reduced, whenever a particular solution of that system is known. A general algorithm for deriving 
the particular solution which makes the reduction possible as been obtained by Xanthopulos |lj]. Holzhey and 
Wilczck found that the Xanthopulos solution one gets from the reduction of the polar equations of Schwarzschild 
and Reissner-Nordstrom black holes are pure gauge solutions i.e. gauge equivalent to the null perturbation. The 
general form of the metric (Q) has, indeed, a gauge degree of freedom, i.e. there exists a one parameter class of 
coordinate transformations which leaves the form of the metric unchanged. In the case of the GHS black hole the 
particular pure gauge solution which reduces the system (|l^)-(|2^) is fl3]| 



X (0 



B. 



(0 



23 
$(0 



-a 2 e-^ 



r + e iv v, r (re-»), r , 
(e "(fe~ v ) tr ) !r — e 4ly i/, r (fe~ ly ) jr , 
e r r r (re ) jr , 



2f 



Following Holzhey and Wilczek, we now introduce a new variable S replacing L 

S = B23 — L — X, 

and make the following substitutions: 



N = N {a) s + n, B 



23 



X = X^s + ^, S = S<®8, 
r 



(24) 
(25) 
(26) 

(27) 

(28) 
(29) 

(30) 
(31) 



The system is now of order six with respect to the new variables n, x, s, b and p since s does appear through its 
derivatives only. It is easy to show from Eq. ( |l9|) that the first derivative of s can be written as a linear combination 
of b, x and p as follows 





FIG. 2. The I = 2 polar potentials and V 2 P (upper panel) and Kf (lower panel) are plotted versus the tortoise coordinate 
for different values of Q = Q e /M. When the charge is set to zero V" reduces to the Zerilli potential, whereas V 2 P and V 3 P are 
the potentials of the wave equation describing pure electromagnetic and scalar perturbations on a Schwarzschild background. 



Similarly, using (pp| ) and (p2|), we can eliminate n in favour of x, b, p and their first derivatives. We are now left 
with only three perturbation variables, x, b and p, governed by three second order equations, namely (pl|), ( |l8| ) and 
(p3|). In (|3l|), x, b and p have been defined in such a way that their first derivative with respect to r* disappears from 
the equations they satisfy. Therefore the system has now the desired form 



drl 



Av, 



where 




(33) 



(34) 



and the components of the symmetric matrix A are complicated functions of r. 
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It is convenient to cast A in the following form: 



A(r) = — [G(r)T(r) + P(r)I], (35) 

where the polinomials D(r), G(r) and P(r) and the matrix T(r) are given in the Appendix, and I is the identity 
matrix. It is remarkable that the eigenvectors of T(r) and A(r) are independent of r. This means that the system 
can be decoupled with a basis transformation in the x-b-p space. The three eigenvalues are the potentials of three 
independent Schrodinger equations: 

d 2 



^ Zf = VfZ? (i = 1, 2, 3). (36) 



It should be noted that for Q e = the equations are already decoupled since the off-diagonal terms of A vanish 
(see Appendix) while, as already pointed out, the GHS solution becomes the Schwarzschild solution. In this case, 
as expected, we find the well-known potentials that rule pure gravitational, electromagnetic, and massless scalar 
perturbations on a Schwarzschild background, respectively. 

Proving the independence of r of the eigenvectors of T(r) directly, even by means of a symbolic procedure, turns 
out to be awkward. A straighforward manner is instead to expand the polinomial valued matrix as: 



T = ^Vt«, (37) 



i=0 



where are matrices independent of r. It is now very easy to show, for instance with the aid of MATHEMATICA, 
that the commute among themselves, so that they can be diagonalized by the same linear transformation. 



The eigenvalue problem for A(r) is now considerably reduced by means of (35) and fl37|). In fact, the three potentials 
can easily be obtained by fixing first the values of the parameters M, Q e and p, then calculating the eigenvalues of 
each Tw, which are now real numbers, and finally summing them up as coefficients of their respective powers of r. 
The eigenvectors of A(r) are then obtained straighforwardly by means of (p5[ ). The three potentials obtained with 
the above procedure appear in a rather simple form as functions of r only. This fact turns out to be extremely useful, 
since to calculate the quasinormal frequencies by means of WKB method one needs up to their sixth derivatives. The 



potentials Vf associated to the wavefunctions Zf (see Eq. 33) are plotted in Fig. || for different values of the 
charge. As in the axial case, when the charge approaches the extremal value, all potentials assume the form of a step. 
When Q e — > 0, all off-diagonal components of the matrix A vanish, and the system of equations ( p3| ) decouple 

into three wave equations for the functions (x,b,p). From eqs. (|3l| ) it is easy to check that, since and 
vanish, the functions b and $ become purely electromagnetic and scalar, respectively. Conversely, the variable 
x is not purely gravitational, since it contains £>23- If, however, we explicitly write the function x in terms of 



B23, X and L in the first Eq. (33), the electromagnetic terms disappear by virtue of the equation satisfied by 6, 



and we are left with the Zerilli equation. 



V. THE QUASI-NORMAL MODES 



We have computed the complex frequencies of the quasi-normal modes associated to the axial and polar wave 
equation by using a WKB approximation devised by Schutz and Will jig] and extended to higher orders by Iyer 
and Will |l6|] . This approach has been applied to the Schwarzchild |l7|] and Reissner Nordstrom |Q cases, and 
for the fundamental quadrupole mode I = 2 agrees with other approaches JTs| ] within 1% both for the real and the 
immaginary parts of the first 3-4 modes. The agreement improves with increasing angular harmonic and decreasing 
mode number. 

The results of our calculations are shown in table | and |J, where we tabulate the real and imaginary part of the 
frequencies of the first five quasi normal modes associated to the axial and polar potentials, for different values of the 
harmonic index £, and of the charge Q e . 

The general behaviour of the quasi-normal frequencies for increasing values of the charge is well described in Figure 
||. There we plot the real and the imaginary part of the first eigenfrequency of the I = 2 mode, associated to the 
potentials V 2 ° and Vf, as a function of the charge Q e . It should be reminded that in the limit Q e — 0, these 
potentials reduce to the Regge- Wheeler and to the Zerilli potentials for the axial and polar gravitational perturbations 
of a Schwarzschild black hole, respectively. It is well known that the axial and polar perturbations of a Schwarzschild 
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black hole is a isospectral, and this is the reason why, in the limit Q e = 0, the axial (V 2 °) and polar (Vf) frequencies 
shown in Fig. ^ converge to a unique value. 

The upper value of Q e /M we consider in our calculations is Q e /M = 1.4, close to the limiting value Q e /M = ^/2, 
where the potential barrier becomes a step which reflects all incident waves. 



TABLE I. The frequencies of the first five quasi normal modes associated to the axial and polar wave equations, are written 
down in units of 1/M, for / = 2 and for increasing values of the charge Q e . 



Vi° 




v 2 a 




yP 






Vi 




vi 




Re(cr) 


Im(cr) 


Re(tr) 


Im(cr) 


Re(a) 




Im(cr) 


Re(cr) 


Im(cr) 


Re(tr) 


Im(cr) 










I = 2 


Qe 


= 0.2 










0.374 


0.089 


0.462 


0.095 


0.374 




0.089 


0.457 


0.095 


0.492 


0.097 


0.347 


0.275 


0.441 


0.292 


0.346 




0.275 


0.436 


0.291 


0.472 


0.297 


0.304 


0.472 


0.408 


0.497 


0.302 




0.471 


0.402 


0.496 


0.442 


0.506 


0.249 


0.674 


0.367 


0.708 


0.244 




0.674 


0.360 


0.706 


0.403 


0.719 


0.180 


0.880 


0.317 


0.921 


0.173 




0.880 


0.309 


0.918 


0.358 


0.935 










I = 2 


Qe 


= 0.4 










0.378 


0.090 


0.479 


0.096 


0.377 




0.090 


0.461 


0.095 


0.516 


0.099 


0.351 


0.276 


0.458 


0.295 


0.350 




0.276 


0.440 


0.292 


0.497 


0.301 


0.308 


0.473 


0.427 


0.502 


0.305 




0.473 


0.406 


0.497 


0.468 


0.512 


0.253 


0.676 


0.387 


0.714 


0.248 




0.676 


0.365 


0.707 


0.432 


0.728 


0.185 


0.883 


0.339 


0.929 


0.177 




0.883 


0.315 


0.920 


0.389 


0.945 










I = 2 


Qe 


= 0.6 










0.386 


0.090 


0.508 


0.098 


0.384 




0.090 


0.474 


0.096 


0.552 


0.101 


0.360 


0.279 


0.489 


0.300 


0.356 




0.278 


0.453 


0.294 


0.535 


0.307 


0.318 


0.477 


0.459 


0.510 


0.312 




0.476 


0.421 


0.501 


0.508 


0.521 


0.264 


0.681 


0.422 


0.725 


0.255 




0.680 


0.381 


0.712 


0.476 


0.740 


0.197 


0.889 


0.378 


0.942 


0.184 




0.889 


0.333 


0.926 


0.437 


0.960 










I = 2 


Qe 


= 0.8 










0.401 


0.092 


0.553 


0.101 


0.396 




0.091 


0.500 


0.098 


0.607 


0.104 


0.376 


0.282 


0.536 


0.308 


0.369 




0.281 


0.481 


0.299 


0.592 


0.315 


0.335 


0.483 


0.509 


0.522 


0.325 




0.482 


0.450 


0.508 


0.568 


0.533 


0.283 


0.689 


0.476 


0.741 


0.269 




0.689 


0.413 


0.722 


0.539 


0.756 


0.220 


0.898 


0.437 


0.962 


0.200 




0.899 


0.368 


0.939 


0.506 


0.980 



I = 2 Qe = 1 



0.429 


0.094 


0.627 


0.105 


0.418 




0.093 


0.545 


0.100 


0.688 


0.107 


0.406 


0.288 


0.612 


0.318 


0.392 




0.286 


0.528 


0.306 


0.675 


0.324 


0.369 


0.492 


0.589 


0.537 


0.350 




0.490 


0.500 


0.520 


0.655 


0.547 


0.322 


0.701 


0.561 


0.762 


0.297 




0.700 


0.466 


0.738 


0.631 


0.775 


0.264 


0.913 


0.529 


0.988 


0.232 




0.913 


0.426 


0.958 


0.603 


1.004 










I = 2 


Qe 


= 1.2 










0.486 


0.097 


0.762 


0.109 


0.466 




0.096 


0.644 


0.105 


0.848 


0.111 


0.468 


0.297 


0.752 


0.330 


0.445 




0.294 


0.631 


0.318 


0.840 


0.336 


0.439 


0.503 


0.736 


0.556 


0.413 




0.502 


0.610 


0.537 


0.826 


0.565 


0.405 


0.715 


0.717 


0.787 


0.375 




0.714 


0.585 


0.761 


0.810 


0.799 


0.365 


0.929 


0.696 


1.019 


0.328 




0.929 


0.556 


0.987 


0.793 


1.034 










I = 2 


Qe 


= 1.4 










0.704 


0.085 


1.252 


0.095 


0.648 




0.083 


1.014 


0.093 


1.403 


0.096 


0.716 


0.258 


1.259 


0.286 


0.666 




0.253 


1.021 


0.280 


1.409 


0.288 


0.744 


0.438 


1.274 


0.480 


0.706 




0.434 


1.038 


0.471 


1.423 


0.483 


0.790 


0.624 


1.297 


0.676 


0.771 




0.624 


1.067 


0.666 


1.444 


0.679 


0.855 


0.816 


1.331 


0.875 


0.863 




0.823 


1.107 


0.865 


1.474 


0.879 







For a Reissner-Nordstrom black hole the perturbed axial and polar equations can be decoupled in terms of four 
functions Z^ and Z^ , respectively ftp] . In the limit Q e = 0, the functions Z^ reduce to the Regge- Wheeler 
and to the Zerilli functions, whereas Z^ reduce to pure axial and polar electromagnetic functions. Moreover, the 
potentials governing Z^ and Z^ are isospectral, as well as those for Z^ and Z^ ■ In Figure |] the same data 
of Figure are compared with the lowest quasi-normal mode frequency of the functions Z^ of a Reissner-Nordstrom 
black hole. Since the extremal value of the charge for a Reissner-Nordstrom black hole is Q e = 1, the data stop at 
that value. 



TABLE II. The frequencies of the first five axial and polar quasi normal modes in units 1/M are tabulated for / = 3 and 
1 = 4, for some values of the charge Q e 



Vi 




v 2 ° 




yv 




vi 




vi 




Re(a) 


Im(o-) 


Re(o-) 


Im(cr) 


Re(<r) 


Im(cr) 


Re(cr) 


Im(cr) 


Re(cr) 


Im(cr) 










I = 3 


Qe = 0.2 










0.601 


0.093 


0.664 


0.096 


0.601 


0.093 


0.657 


0.096 


0.687 


0.097 


0.584 


0.282 


0.649 


0.291 


0.584 


0.282 


0.642 


0.290 


0.672 


0.294 


0.555 


0.477 


0.623 


0.492 


0.555 


0.477 


0.616 


0.490 


0.647 


0.497 


0.517 


0.678 


0.590 


0.698 


0.518 


0.678 


0.582 


0.696 


0.615 


0.704 


0.473 


0.882 


0.550 


0.907 


0.473 


0.882 


0.542 


0.904 


0.578 


0.915 










t = 3 


Qe = .4 










0.607 


0.093 


0.686 


0.097 


0.606 


0.093 


0.666 


0.096 


0.716 


0.098 


0.590 


0.283 


0.672 


0.294 


0.589 


0.282 


0.651 


0.291 


0.702 


0.298 


0.562 


0.479 


0.647 


0.496 


0.560 


0.478 


0.625 


0.492 


0.678 


0.503 


0.525 


0.680 


0.614 


0.704 


0.523 


0.680 


0.592 


0.698 


0.648 


0.712 


0.480 


0.885 


0.577 


0.915 


0.479 


0.885 


0.552 


0.908 


0.612 


0.926 










I = 3 


Qe = .6 










0.621 


0.094 


0.724 


0.099 


0.617 


0.094 


0.688 


0.097 


0.762 


0.100 


0.605 


0.285 


0.711 


0.299 


0.601 


0.284 


0.673 


0.294 


0.749 


0.303 


0.576 


0.482 


0.687 


0.504 


0.572 


0.481 


0.648 


0.497 


0.726 


0.511 


0.540 


0.685 


0.657 


0.715 


0.536 


0.684 


0.616 


0.705 


0.698 


0.724 


0.497 


0.891 


0.622 


0.928 


0.492 


0.889 


0.579 


0.916 


0.665 


0.940 










I = 4 


Qe = 0.2 










0.811 


0.094 


0.862 


0.096 


0.811 


0.094 


0.854 


0.096 


0.882 


0.097 


0.799 


0.285 


0.851 


0.290 


0.798 


0.285 


0.842 


0.289 


0.870 


0.292 


0.776 


0.479 


0.801 


0.691 


0.776 


0.479 


0.821 


0.487 


0.849 


0.492 


0.746 


0.679 


0.768 


0.898 


0.745 


0.679 


0.792 


0.690 


0.822 


0.696 


0.710 


0.882 


0.729 


1.107 


0.709 


0.882 


0.758 


0.895 


0.789 


0.904 



I = 4 Qe = 0.4 



0.820 


0.095 


0.890 


0.097 


0.818 


0.094 


0.868 


0.096 


0.916 


0.098 


0.808 


0.286 


0.879 


0.293 


0.806 


0.285 


0.856 


0.291 


0.905 


0.296 


0.785 


0.481 


0.858 


0.493 


0.783 


0.481 


0.835 


0.490 


0.885 


0.498 


0.755 


0.681 


0.831 


0.698 


0.753 


0.681 


0.807 


0.693 


0.858 


0.704 


0.720 


0.885 


0.798 


0.906 


0.718 


0.884 


0.774 


0.899 


0.827 


0.913 










I = 4 


Qe = 0.6 










0.841 


0.095 


0.936 


0.099 


0.835 


0.095 


0.899 


0.097 


0.969 


0.100 


0.829 


0.288 


0.926 


0.298 


0.823 


0.287 


0.888 


0.294 


0.959 


0.301 


0.807 


0.485 


0.906 


0.501 


0.801 


0.484 


0.867 


0.495 


0.941 


0.506 


0.778 


0.686 


0.881 


0.708 


0.771 


0.685 


0.840 


0.700 


0.916 


0.715 


0.743 


0.891 


0.850 


0.919 


0.736 


0.890 


0.808 


0.908 


0.887 


0.927 
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Re(a) lm(cr) 




FIG. 3. The real and imaginary part of the frequencies of the quasi-normal modes of the dilaton black holes are plotted 
versus the electric charge Q — Q e /M, for £ = 2, and for the two potentials which reduce to the Regge- Wheeler and to the 
Zerilli potentials when the charge is set to zero, i.e. V-f (triangles) and Vf, (circles). 

We see that the real part of the frequency increases as a function of the charge both for a dilaton, and for a 
Reissner-Nordstrom black hole; similarly, the imaginary part increases and then decreases to a finite value as the 
charge approaches the limiting value. Since however since the limiting values are different, the imaginary part for 
Reissner-Nordstrom begins to decrease while the corresponding GHS values still increase. 

Figures || and || clearly show that the potentials V 2 ° and Vf of a dilaton black hole are not isospectral as they 
are instead when Q e = 0. In addition, it should be reminded that the potential V$ rules the equation for Z% 
which is a combination of electromagnetic and gravitational perturbations only (cfr. Eq. |i~3| ), whereas Vf appears 
in the polar equation for Zf, which is a combination of electromagnetic, gravitational and scalar perturbations. 



Re(cr) lm(cr) 




FIG. 4. The eigenfrequencies of the quasi-normal modes of the dilaton black hole given in Figure 3, are compared with 
those of a Reissner-Nordstrom black hole for I = 2. For Q e = both solutions reduce to Schwarzschild and the lines converge to 
its first, I = 2, pure gravitational mode a — 0.373 + 0.089i. For increasing values of the charge Reissner-Nordstrom frequencies 
remain isospectral, while GHS frequencies split in two parts, depending on whether they belong to axial or polar modes (see 
text). Since the extremal value of the charge for a Reissner-Nordstrom black hole is Q e /M = 1, the data stop at that value. 
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Finally, in Figure |g we plot the real and imaginary part of the frequency of the lowest I — 2 quasi-normal mode for 
the remaining potentials, V", Vl? and V£ ', as functions of the charge Q e . When Q e = 0, V\ a and V 2 P reduce 
to those governing the pure electomagnetic, axial and polar perturbations of a Schwarzschild black hole, respectively, 
and indeed, they are isospectral in that limit. In the same limit, V£ reduces to the potential of the pure scalar 
perturbations of a Schwarzschild black hole. 



Re(cr) lm(cr) 




FIG. 5. The real and imaginary part of the frequencies of the lowest I = 2-quasi-normal modes of the dilaton black hole, are 
plotted versus the electric charge Q — Q e /M, for the three potentials V\ (triangles), V% (circles) and V 3 P (stars). In the limit 
Q = 0, Vi and V 2 P reduce to the potentials governing the axial and polar pure electromagnetic perturbations of a Schwarzschild 
black hole, whereas VJf becomes the potential which appears in the wave equation of the pure scalar perturbations. 

VI. CONCLUDING REMARKS 

The spectrum of the quasi-normal modes of a charged, dilaton black hole is different from that of a Schwarzschild 
or a Reissner-Nordstrom black hole. For a Schwarzschild black hole the perturbations are completely described by 
the Rcggc- Wheeler and the Zcrilli equations for the axial and the polar perturbations. Although the analytic form 
of the two potential barriers is different, they are related by a very simple equation (MT, ch. 5, & 43) which allows 
them to have the same reflexion and trasmission coefficients. Therefore, since the quasi-normal mode frequencies are 
the singularities of the scattering amplitude, it follows that the two potentials are isospectral. Thus, a perturbed 
Scharzschild black hole emits axial and polar gravitational waves at exactly the same frequencies. 

For a Reissner-Nordstrom black hole the perturbed equations can be reduced to four wave equations, two for the 
axial and two for the polar perturbations, respectively. The four wave- functions Zf 1 and Z 2 , where + stands for 
polar and — for axial, are a linear combination of the gravitational and the electromagnetic functions belonging to the 
corresponding parity. It turns out that the two potentials and Vy are again related in such a way that they 

have the same reflection and trasmission coefficients, and the same is true for VJj . Thus the coupling GW-EM is 
such that it preseves the isospectrality of the axial and polar perturbations. However there is an important difference 
with respect to the Schwarzschild case: no quasi-normal mode exists that is purely electromagnetic or gravitational, 
which means that the excitation of a mode will be accompanied by the simultaneous emission of both gravitational 
and electromagnetic waves. 

For a charged black hole in a theory described by action ([j]) the situation is different. Let us consider the axial 
perturbations first. As shown in Section 3, the perturbed equations can be reduced to two wave equations, but 
the perturbed dilaton does not couple to the electromagnetic and gravitational fields. This is due to the fact that 
the dilaton is a scalar, and its axial perturbation vanishes. Consequently, the excitation of an axial mode will be 
accompanied only by the emission of gravitational and electromagnetic waves. However, the dilaton appears in the 
unperturbed metric functions that determine the shape of the potentials of the axial wave equations. Thus it affects 
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the scattering properties of the axial potentials. It is interesting to note that the real part of the quasi-normal mode 
frequencies of the axial potential are very similar to those of the Reissner- Nordstrom black hole (see Fig. ||), 

even though the dilaton solution does not reduce to the Reisnner-Nordstrom solution in any limiting case. 

On the other hand, the two wave equations which describe the polar perturbations of a Reissner-Nordstrom black 
hole, are replaced by three wave equations in the case of a dilaton black hole, and they couple the gravitational, elec- 
tromagnetic and scalar perturbations. This occurrence breaks the symmetry between axial and polar perturbations, 
and makes the scattering properties of the two parities different (see Fig. ||) . 

We conclude that for a dilaton black hole the excitation of an axial mode induces the simultaneous emission 
of gravitational and electromagnetic waves, whereas the excitation of a polar mode is accompanied by the further 
emission of scalar radiation. In addition, gravitational and electromagnetic polar waves are emitted with frequencies 
and damping times different from the axial ones. 



ACKNOWLEDGMENTS 



F.P. whishes to thank Roberto De Pietri, Antonio Scotti and Michele Vallisneri for useful discussions and help in 
computing technicalities. 



VII. APPENDIX 



In this Appendix we provide all the elements of the matrix A that appears in equation 
contained all the relevant features of the polar perturbations. It proves convenient to cast A in t 



A(r) 



1 



D(r) 



[G(r)T(r) + P(r)I], 



and in which are 
ic following form 

(38) 



where I is the identity matrix, and D(r), G(r) and T(r) are: 



D(r) 



4pQ e r 4 
x(3M( 
x(4MQ e 4 - 14 Af 2 Q 



Qe 2 + Mr) 
- 7M 2 Q e 2 r 



Q e 4 r + 6 Af 3 r 2 - Af p 



2 r 2 



M 2 p 2 r 3 ) 



Q 4 r + 12 Af r-2M/i Q 2 r +2 Af p 2 r 3 ) 



G(r) = pQ e (2M - r) r 2 x (-Q e 2 + Mr) 
(3MQ e 4 - 

P{r) = (1936 M e Q 



7M 2 Q 2 r-Q 4 r 



6 M 3 r 2 — Af [i 2 Q 2 r 2 + M 2 p 2 r 3 ) /Af, 



8 



2432 Af 7 Q 4 r - 448 M b Q e r r 



324 M Q e 

f 5n 6 



160 M 4 n 2 Q e 8 + 
f 624 Af 5 fi 2 Q e 6 ', 



:M 2 Q, 



10 



Q. 



12 



2112 Af 8 Q 2 r 2 



1136 M 6 Q e 4 r 2 



1520 AT 6 p? Q e 4 r 2 



\M A n z Q e a r + 
64M 4 Q e 6 r 2 



232 Af 4 p 2 Q e ° r z + 24 AT 4 ^ Q e ° r z - 8 Q 



6 Jl 



f4 ,4, 



6 J2, 



f2 ,,2 , 



8 „2 



1152 A/ 9 r 3 + 2688 AT 7 Q 2 r 3 



1984 M 7 n 2 Q 2 r 3 



4„3 



272Af D /x Q e 4 r' 
-4 AT 3 (j 4 Q e 6 r 3 - 1152 AT 8 r 4 



r 

2 „4 



2 r 3 

1 16 A/ 3 /' 2 ' 

r8 ,,2 A 



264 Af ^ Q e 4 r 3 + 16 Af p z Q, 



6 r 3 



464 AfV Q e ~r 
16Af 4 M 6 Q e 4 r 4 
32 Af 5 fi 6 Q 



2 r 5 



960 Af 8 n A r 4 + 832 Af b Q t 
- 16 M 4 v 2 Q 4 r 4 - 28 Af 4 fi 4 Q 4 r 4 - 
384 Af 7 [i 2 r 5 - 224 Af 7 p 4 r 5 + 64 M 5 fi 4 Q c 
32 Af 6 [i 4 r 6 - 16 AfV 6 r 6 )/8M 2 fiQ e . 



2 r 5 



Tn = 2448 Af b Q e D + 324 A/ 4 Q e s - 160 Af 4 p 2 Q 



iAf 2 Q e 



Qe 



8Af 3 Q e 



(-496 Af 4 - 64 Af 2 Q 2 + 110 Af 2 pt 2 Q 2 - 4 Q 4 + 11 ^i 2 Q e 4 ) r + 8M 2 > 
(408 Af 6 + 42 Af 4 Q e 2 - 222 Af 4 p 2 Q 2 - 45 Af 2 p 2 Q 4 + 3 Af 2 p 4 Q 4 - 
/i 2 Q e 6 ) r 2 + 4 Af 3 (-288 Af 6 + 400 Af 4 p 2 Q 2 + 108 Af 2 p 2 Q 4 - 
34 Af 2 p 4 Q 4 + Ap 2 Q e e - p 4 Q e 6 ) r 3 + 4 Af 4 ^ 2 (-144 Af 4 - 
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32 M 2 Q 2 + 52 M 2 y 2 Q 2 - 4 Q 4 - 7 (J 2 Q 4 - 4 ^ 4 Q e 4 ) r 4 + 
32 M 5 M 4 (-3 M 2 + 2 Q e 2 + y 2 Q e 2 ) r 5 - 16 M 6 y 4 (2 + ^ 2 ) r 6 

T 2 i = 128 M 5 y Q e 7 + 8 M 2 y Q e 5 (-76 M 4 + 8 M 2 y 2 Q 2 - 

Q e 4 ) r + 16 M 3 n Q 3 (56 M 4 + 22 M 2 Q 2 — 8 M 2 y 2 Q 2 + Q 4 - 

2 M 2 Qe 4 ) r 2 + 32 M 4 /z Q e (-12 M 4 - 24 M 2 Q 2 + 

2M 2 M 2 Q e 2 + 6^ 2 Q e 4 + M 4 Q e 4 )r 3 + 32M f >Qe 

(12 M 2 - 9 m 2 Qe 2 - 2 M 4 Q e 2 ) r 4 + 32 M 6 M 3 (4 + y 2 )Q e r 5 

T 3 i = 16 M 3 it V2 + /i 2 Qe 2 r (Q e 2 - Af r) (14 M 2 g e 2 - Q 4 - 

24 M 3 r + 4I/i 2 Q e 2 r + 6 M 2 r 2 - 6 M 2 y 2 r 2 - ^ 2 Q e 2 r 2 +2M/i 2 r 3 ) 

T 22 = 20 64 M 6 Q e 6 + 260 M 4 Q e 8 - 160 M 4 ^ 2 Q e 8 + 8 M 2 Q e 10 + Q e 12 + 
4MQ e 4 (-536 M 6 - 164M 4 Q e 2 + 116 M 4 y 2 Q 2 + 
6 M 2 Q 4 + 18 M 2 y 2 Q 4 + Qe 6 ) r + 8 M 2 Q 2 
(72 M 6 + 6 M 4 Q 2 - 98 M 4 y 2 Q 2 - 32 M 2 Q 4 - 

43 M 2 /i 2 Q e 4 -5M 2 ,! 4 Qe 4 - Qe 6 + ^ Qe 6 ) r 2 + 

4 M 3 Q e 2 (432 M 4 + 256 M 4 yu 2 + 48 M 2 Q 2 + 

12 M 2 y 2 Q 2 - 18 M 2 (i 4 Q 2 - 12 y 2 Q e 4 - 

M 4 Qe 4 ) r 3 +A M 4 (-288 M 4 - 144 M 4 y 2 + 160 M 2 y 2 Q 2 + 

76 M 2 y 4 Q 2 + 8 y 2 Qe 4 - 11 M 4 Qe 4 - 

4 M 6 Qe 4 ) r 4 + 16 M 5 y 2 (-24 M 2 - 12 M 2 y 2 + 5 ^ 2 Q 2 + 

2 M 4 Qe 2 ) r 5 - 16 M 6 jU 4 (2 + M 2 ) r 6 



T 32 = 8 M 2 v / 2~+7t 2 Qe (-16 M 3 Qe 6 - 36 M 4 Q 4 r + 

8M 2 Q e 6 r~ 8 M 2 ^i 2 Q e 6 r + Q e 8 r + 192 M 5 Q e 2 r 2 - 

52 M 3 Q e 4 r 2 - 16 M 3 ^ 2 Q e 4 r 2 - 2 M Q e 6 r 2 + 

4M/i 2 Q e 6 r 2 - 144M 6 r 3 +48 M 4 Q 2 r 3 + 72M 4 fJ 2 Q 2 r 3 - 

16 M 2 y 2 Q 4 r 3 — 4M 2 M 4 Q e 4 r 3 - 48 M 5 M 2 r 4 + 

12 M 3 M 2 Q 2 r 4 + 8M 3 y 4 Q e 2 r 4 — 4 M 4 M 4 r 5 ) 

T 33 = 1936 M 6 Q e 6 + 324 M 4 Q e 8 - 160 M 4 y 2 Q e 8 + 8 M 2 Q e 10 + Q 42 + 
8 M 3 Q 4 (-304 M 4 - 56 M 2 Q 2 + 78 M 2 y 2 Q 2 + 
11 y 2 Q e 4 ) r + 8 M 2 Q 2 (264 M 6 - 142 M 4 Q 2 - 
190 M 4 [i 2 Q 2 — 8 M 2 Q 4 - 29 M 2 Q 4 + 

3 M 2 (i 4 Q 4 - M 2 Qe 6 ) r 2 + 4 M 3 (-288 M 6 + 672 M 4 Q e 2 + 
496 M 4 y 2 Q 2 - 68 M 2 fj 2 Q 4 - 66 M 2 fi 4 Q 4 + 

4 M 2 Q e 6 - y 4 Q e 6 ) r 3 +4M 4 (-288 M 4 - 240 M 4 y 2 + 
208 M 2 y 2 Q 2 + 116 M 2 y 4 Q 2 -Ay 2 Q 4 - 

7 m 4 Qe 4 - 4 ^ 6 Q e 4 ) r 4 + 32 M 5 /i 2 (—12 M 2 — 

7 M 2 y 2 + 2 y 2 Q 2 + M 4 Q e 2 ) r 5 - 16 M 6 M 4 (2 + M 2 ) r 6 
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